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, The long-range energy deposition by heavy-ion beams makes new types of shock 

I wave experiments possible in the laboratory. We have investigated a situation that 

is of relevance to supernova dynamics in astrophysics, where a shock wave is ir- 
radiated by a flux of neutrinos depositing energy throughout the shock wave and 
surrounding matter, thus changing the behaviour of the running shock. We have 
' carried out fluid-dynamical simulations to study generic features of stimulated shock 

. waves. First we consider an idealized case assuming uniform energy deposition into 

I a planar shock wave propagating through an ideal gas. Then we investigate more 

realistic situations realizable in laboratory experiments with heavy-ion beams. We 
have found that energy deposition leads to two important effects: acceleration of 
the shock front and decay of the shock strength. The possibility of laboratory 
experiments is briefly discussed. 
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?^ 1 Introduction 



Supernova explosions represent one of the most spectacular phenomena in our Universe. 
In recent years much effort has been devoted to developing realistic theoretical models 
of this complex process (see the recent review [|l|). The main mechanism includes the 
shock wave generation upon the bounce of the infalling iron core of a massive star, but 
detailed calculations showed that the prompt bounce shock mechanism does not lead 
directly to ejection of the stellar envelope. Due to severe energy losses due to, e.g., 
photodisintegration of iron nuclei, the outward shock wave stops inside the iron core. 
Recently several new mechanisms of the shock revival have been suggested 0, |^ which 
include neutrino heating and convection in the postshock matter. These processes raise 
the postshock pressure and provide additional energy for the shock wave expansion. The 
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calculations performed by various groups differ, however, in conclusions concerning the 
possibility of successful shocks. In this rather unclear situation a better understanding 
of the stimulated shock dynamics is highly desirable 0. Fortunately, intense heavy-ion 
beams open a new possibility to study such processes in laboratory experiments. 

In the first part of this paper we formulate a simple model to study the evolution 
of a planar shock wave under influence of homogeneous irradiation. We have performed 
detailed fluid-dynamical calculations for the ideal gas equation of state. The matter flow 
behind the shock front is rather complicated, resembling a decay of an initial disconti- 
nuity 0]. The calculations show that under irradiation the shock front accelerates, but 
the density jump at this front diminishes. It is remarkable that a self-similar regime of 
hydrodynamic flow is established at late times. 

In the second part of the paper we generalize our model to consider more realistic 
situations which may be easily reproduced in laboratory conditions. Namely, we consider 
the case when the initial shock wave is created by a beam with finite extension in the 
transverse directions. Specifically, we study the scenario in which the first beam creates 
a shock wave in a cylindrical target and then, when it has developed sufficiently, the 
target is irradiated by a second beam. The target length is chosen in such a way that 
the shock generated at the initial stage is completely within the deposition region of 
the second beam[|. In our simulations we assume that both are heavy-ion beams with 
equal properties and are directed along the target axis from opposite directions. General 
properties of hydrodynamical flow in such targets have been studied earlier in Ref. ||^. 
Another interesting phenomenon is predicted for the case of a single beam irradiation 
with constant temporal profile. It is shown that with increasing irradiation time, when 
radial flow of matter behind the shock front becomes noticeable, the shock wave in its 
central parts decays as compared to initial stages. The origin of this phenomenon consists 

in the increased range of the bombarding particles at later stages of target irradiation. 

"'^It is, of course, irrelevant how the shock itself is created — instead of a heavy-ion beam a laser could 
be used, for example, — and the second beam could come from any suitable direction. 
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The article is organized as follows: in Sect. II a simple fluid-dynamical model is 
formulated for the case of homogeneous energy deposition. Then in Sect. Ill this model is 
used to study the dynamics of a planar shock wave in matter with an ideal gas equation of 
state. The asymptotic regime of the matter flow is considered by a semi-analytic method 
outlined in the Appendix. The results of more realistic calculations and suggestions for 
future experiments with heavy-ion beams are given in Sect. IV. The main results of the 
present paper are summarized in Sect. V. 

2 Equations of fluid dynamics 

Let us consider a beam of energetic particles irradiating a target with the mass density p . 
Often we shall speak about photon irradiation, but the same approach can be applied 
for other penetrating particles, e.g. for neutrino beams. The strength of irradiation is 
characterized by the energy deposited per unit space-time volume in the local rest frame 
of the target, 

where A is the specific deposition rate. Let the target consist of atoms of different kinds i 
with particle number densities rii . In the case of monochromatic photons with the energy 
flux /(a;) one has 

Xp^I{u;)j:n,<af\u;)>, (2) 

i 

where hcu is the photon energy, < af'^ > is the cross section of photon absorption on 
atoms of the i-th kind. The angular brackets denote averaging over their momentum 
distribution. In this paper we neglect internal heat transport processes and assume A to 
be constant in time and space. The last assumption implies that the external radiation 
is homogeneous, Doppler effects are not important and the photon absorption length, 
{J2ni < af'^ >)^) is large as compared to the characteristic size of a target. 

i 

Assuming that energy deposition proceeds under condition of local thermodynamic 
equilibrium, we study the dynamics of target flow by solving nonrelativistic fluid-dynamical 



equations. Let us consider first a one-dimensional case when matter moves along the x 
axis with the velocity v = Vx ■ The equations of motion then take the form 

dtp + dx{pv) = 0, (3) 
dtipv) + dx [P + pv^) = , (4) 
dt (e + pvy2) +dx[v[e + P + pt^'/s)] = Ap . (5) 

Here P is pressure and e is the energy density. By using Eqs. and thermodynamic 

relations one may rewrite Eq. (^) in the equivalent form: 

T{dt + vdx)sm = X, (6) 

where T is temperature and Sm is the specific entropy. This equation shows that A is, in 
fact, the specific rate of heat deposition into the target matter. 

In the following we assume that the target matter can be regarded as an ideal Boltz- 
mann gas with one species of molecules and a constant ratio of heat capacities 7 = Cp/C^. 
The equation of state of such matter has a very simple form 

P={^-l)e = pRT/M, (7) 

where R is the gas constant and M is the molar mass. The specific entropy in this case is 

R P 

s^ = j--^\n- + const. (8) 

Using Eq. (|^ and introducing the specific volume V = 1/p one can rewrite Eqs. 
as follows 

{dt + vdx)V = VdxV, (9) 
{dt + vdx)v = -VdxP, (10) 
V{dt + vdx)P + iPidt + vdx)V = {^-1)X. (11) 

It is convenient to introduce instead of P, V, v new dimensionless quantities P = P/Pq, 
V = poV, V = v/vo where Pq and po are initial pressure and density, and 



vo = JPo/po. (12) 



Let us make the transition from t,x to dimensionless variables t = t/to and x = x/xq , 
where 

to = Vq/X, Xo = voto. (13) 

Now Eqs. (^-(^) may be rewritten in the same form, but with the replacements 
P, V, v,t,x P, V, V, i, X and A ^ 1 . If the initial conditions do not impose additional 
scales, there is no need in solving fluid-dynamical equations for different deposition 
rates A . It is sufficient to find the solution only for one flxed value of A , then the above 
scaling can be used to obtain the solution for other values. For example, this scaling can 
be applied for the shock wave initial conditions (see the next section). 

We close the general discussion by defining three families of characteristics 
They are solutions of the following differential equations 

x{t) = v + Cs (C+), (14) 

x{t) = v-cs (C_), (15) 

x{t) = V (Co), (16) 

where v and are taken at x = x{t) . As well-known these characteristics describe 
propagation of small disturbances of fluid-dynamical quantities. In particular, the entropy 
disturbances propagate along the Co characteristics (the latter are also the collective flow 
trajectories). 

In the following we assume that irradiation starts at t = and the initial proflles 
p(x, 0) , f (x, 0) and P(x, 0) are known function of x. Let us consider flrst the case of 
homogeneous initial conditions, when p(x, 0) = po; v{x^ 0) = vo and P(x, 0) = Pq. In this 
case the solution of Eqs. (|9|)-(p!TD is trivial: 

V{x,t) = 1/po, (17) 
v{x,t) = Vo, (18) 
P(x,t) = Po + (7-l)PoAt. (19) 
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In the same case the adiabatic sound velocity is equal to 



Cs = JiPV=Jcl + ^{^-l)\t, (20) 



where cq = y'jPo/po is the initial sound velocity. 

For homogeneous initial conditions, substituting ([T8|) , (|20|) into Eqs. (|l4D -(p!^), one 
can find characteristics analytically: 
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x^Jt) = VQt±^-^-—jj^[c'o + ^{^-l)Xt\ + const, (21) 
x^^(t) = vot + const . (22) 

According to Eq. (^), at large irradiation times the C± characteristics approach asymp- 
totically the same lines which do not depend on initial conditions: 

x^Jt)^±eoV^ (t»c^/A,v2/A), (23) 

where 

^0 = ^7(1-1). (24) 

3 One— dimensional shock waves 
3.1 Numerical results 

Let us consider a one-dimensional shock wave propagating to the left along the x axis 
(directed to the right) and denote the position of its front by a;sh(^) • For brevity we 
introduce the vector y = (V, v, P) which combines the set of fiuid-dynamic quantities in 
a compact way. It is assumed that at x = Xgh this vector jumps from yi = y{xsh — 0,t) to 
y2 = y{xsh + 0,t). In the rest-frame of the shock wave (moving with the velocity D = \xsh\ 
with respect to the laboratory frame) the fiuxes of mass, energy, and momentum should 
be continuous at x = Xsh ■ This gives three relations H] connecting D and the components 



6 



of vectors yi^2 ■ In the case of an ideal gas one obtains^ 



vi + D ^ V2 + D ^ P2 -Pi ^ h (Pi + P2) , . 

Vi V2 \V1-V2 V1 + V2 ' ^ ' 

where fi 2 are defined in the laboratory frame. Therefore, at given t/i, all characteristics 

of the shock wave can be determined if one of the quantities P2, V2, V2 is known. One 

should bear in mind that the relations ( pSf ) are based only on the local conservation laws 

and hold also for shock waves with energy deposition. However, as will be shown below, 

in this situation the characteristics of the shock wave (in particular, D) are in general 

t ime-dep endent . 

It is further assumed that a step-like one-dimensional shock wave was created in a 
target at t < . If pressure and density before and behind the shock front do not change 
with X and t, the shock velocity D will be constant. Below we study the dynamics of 
this "initial" shock wave after switching on the energy deposition at t = . Eqs. (^-(|ll|) 
are solved numerically by using the FCT-algorithm 0. The calculations are performed 
in the rest frame of the initial shock. Choosing the origin of the x axis at the position 
of the shock front, one has Xsh{t) = at t < . To be more specific, the following initial 
conditions are applied 

yix,0)=yioe{-x)+y2oe{x), (26) 

where 6(x) = ^(1 + signx) , and yio and 7/20 denote the fluid-dynamical quantities in 
front and behind of the initial shock front, respectively. We choose these quantities by 
fixing the density (po = Pio = l/^io) and pressure (Pq = Pio) at x < as well as the 
pressure (P20) at x > . 

Using the relations (531) one has 



_ (7-l)x + 7 + l 
(7 + l)x + 7- 1 

^10 = V20T^ = D{0) = Vo\ , (28) 

L20 V ^ 



•^It is worth noting that the last equahty in Eq. ( p5| ) is in fact the Rankine-Hugoniot rela- 
tion. Excluding D from the first two relations one can find that the velocity jump is equal to 



Vl~V2^ ^/{P2~ Pl)iVl-V2). 



where x = P20 / Pio is the initial pressure ratio, vq is defined by Eq. (0) and the velocities 
Vio {i = 1, 2) are taken in the rest frame of the shock front. 

Most calculations have been performed for 7 = 5/3 and x = 4 . In this case the initial 
density jumps at the shock front from po to 2.125 po and the fiow velocity jumps from 
fio — 2.380^0 to V20 — 1.120^0- Some results of numerical calculations for these initial 
conditions are shown in Figs. 1-4. The results are presented using the scaled variables 
t/to and x/xq , with to, defined in Eq. (|1^). As explained in the preceding section, in 
such a representation the whole A dependence is contained only in the scales to and Xq. 

Fig. |1] shows the density profiles at different times from the beginning of irradiation. 
One can see that irradiation deforms the shape of the initial shock wave at t > . The 
region of the perturbed flow is wider at later times. With increasing t the shock wave 
becomes weaker (the density jump at the left hand side diminishes). However, the absolute 
value of the shock front velocity becomes larger as compared to the initial value -D(O) . 

The acceleration of the shock front is clearly visible in Fig. |l|. Indeed, the shift of the 
shock front position from the point a; = increases nonlinearly with t . The reason for 
acceleration is rather simple: as follows from Eq. pO|) and Fig. |^ the sound velocity Cg 
behind the shock front increases with t due to the heat deposition. 

At given t the density proflle has three characteristic points: 1) the shock wave discon- 
tinuity at X = Xsh{t) , 2) the intermediate kink where the derivative d^p jumps and 3) the 
boundary of the perturbed region on the right hand side. The position of the third point 
can be found analytically: it is in fact the characteristic x^_^{t) with the initial condition 
x^^(O) = . The calculation conflrms that this point moves in accordance with Eq. ( PT| ) 
where vo = V20 and cq = C20 = 1/7^20^20 • It can be shown that the second point corre- 
sponds to the characteristic x^^(t) with x^^{0) = . The direct numerical integration of 
Eq. (|l^) gives the result shown by the dotted curve in Fig. |l|. One can see that on the 
X — p plane the Cq characteristic indeed goes through the kink positions following from 
Eqs. 
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Profiles of pressure calculated for the same initial conditions are shown in Fig. ^ One 
can see that the pressure jump at the shock front, P2 — Pi, becomes larger at later stages 
of irradiation. We have checked that the density and pressure jumps at a; = x^h satisfy 
the Hugoniot adiabate given by Eq. (^51). In accordance with Eq. ([T9|) , pressure increases 
linearly with t a.t x ^ ±00 . At known p and P the sound velocity Cg can be calculated by 
using the first equality in Eq. (^). In the case of the ideal gas Cg is proportional to VT . 
The results of the calculation are shown in Fig. ^ It is seen that the jump in c^, and 
therefore also in T , practically does not increase with t . At large t the asymptotic values 
of Cs at X ^ ±00 become nearly equal. This agrees with Eq. (pO]). Indeed, at t ^ c|o/A 
one obtains the relations 

lim c, ~^eoV^, (29) 

where ^0 is defined by Eq. (pi]). 

Fig. ^ shows the results for the flow velocity profiles. Similar to density, values of v 
do not change with t at x — ±00 . However, the flow velocities in the perturbed region 
diminish as compared to initial values. As will be shown in the next section these velocities 
become negative at late stages of irradiation. Numerical values of the velocity jump agree 
with conditions given in Eq. (see footnote on page 7). At large t this jump increases 
approximately proportional to \/i . 

3.2 Self— similarity of flow at large irradiation time 

The results presented in Sect. IIIA suggest that at large irradiation time, the fluid dy- 
namical quantities exhibit some properties of self-similarity. At t ^ c^q/ A characteristic 
values of pressure and sound velocity increase, respectively, like t and v^. On the other 
hand, at late stages of irradiation the size of the region with a perturbed flow grows 
approximately like t^^^ . This means that at t ^ 00 the following relations should hold 
asymptotically 

V ^ Ks = Po'V^(0, (30) 
9 



p 



V 



fas = (0 , 

Pas = (7-l)PoAtP(0, 



(31) 



(32) 



where scaled dimensionless quantities V, v, P depends on x and t only via the self-similar 
variable 



By substituting asymptotic expressions (^0|)-(^) into the equations of fluid dynamics 
(^-(0) one can obtain the set of ordinary differential equations for V, v and P (see 
Appendix). Apparently, these equations can not be solved analytically; it is possible, 
however, to find the limiting behaviour of flow directly from the numerical solution of 
Eqs. (D-(0) at large t . 

Fig. 1^ represents density as function of ^ at several fixed times t . One can see that in 
agreement with Eq. (pO]), at large t the density depends on x and t only via ^ . The shock 
wave does not disappear a.st —>■ oo. At late stages its front accelerates in accordance with 
the relation 



where ^sh is a constant fully determining the asymptotic properties of the shock wave. In 
particular, ^sh defines asymptotic behaviour of the shock front velocity: 



As seen in Figs. 5-6, the region of the perturbed (nontrivial) flow corresponds to the 
interval ^sh < ^ < • It will be shown below that the asymptotic position of the shock 
front ^ = ^sh moves to the left with rising initial pressure ratio x ■ For weak initial shocks, 
X ~ 1 , the shock front becomes close to the limiting position of the C_ characteristics, 
i.e. ^sh ^ -^0- 

Another feature of the self-similar regime, is the appearance of a rapid density varia- 
tion at some intermediate behind the shock front. By analyzing Eqs. (|A1|) -(|A^) (see 




(33) 




(34) 




(35) 
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Appendix) one can show that at large t the disturbed zone is characterized by the combi- 
nation of a strong discontinuity (the shock front at ^ = ^sh) and two weak discontinuities^ 
at = and ,^0 • According to Fig. ^ the convergence to the self-similar behaviour is 
much slower for the velocity profiles. It is interesting to note that the flow velocity changes 
sign at late stages of irradiation. Comparing Figs. || and ||, one can see that the point 
^ = ^* corresponds to the minimum of the asymptotic velocity profile v{C,) . This point 
also gives the limiting position of the Cq characteristics at t 00 . Indeed, substituting 
(pl|) into Eq. (|16|), one can see that all Cq characteristics satisfy the asymptotic relation 
x^^{t) ~ ^^VXt^ , where ^* is determined from the equation^ 

v{^*) = l^.. (36) 

Figs. 7-9 show asymptotic profiles of flow velocity, density and pressure calculated for 
initial ratios x = 2, 4 and 10 . In accordance with the above discussion, the line 3^/2 goes 
through minima of the asymptotic velocity profiles at different x ■ This is clearly seen in 
Fig. ^ Asymptotic density profiles are shown in Fig. Here the points C,* are marked by 
arrows. Fig. ^ shows the asymptotic pressure profiles which, unlike the density profiles, 
vary smoothly at ^ = . 

Jumps of asymptotic parameters at the shock front ^ = ^sh can be determined by using 
equations (^) in the limit t 00. At ,^ < ,^sh one has ■0 = and P = V = 1 . Marking 
the flow parameters at ,^ = ,^sh + by indices 'sh', one gets the relations 

3 



^sh - - 6 



sh 



/ , N -Psh - 1 3 
(7 - 1) TT = n^O 



1 + Psh ^3^^ 



Using these conditions, Vsh , Psh and V^h may be obtained as functions of ^sh • For example, 
the scaled pressure behind the shock front is equal to 



7 + 1 



-sh=l + 7^(^-l . (38) 



^In such discontinu ities spatial derivatives of fluid-dynamical variables are infinite or exhibit jumps. 
^According to Eq. (Al), ii is minimal at ^ = . 
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Table 1: Parameters of asymptotic flow behind the shock wave (7 = 5/3) 



X 


^sh 










2 


-0.720 


1.038 


-0.040 


1.062 


-0.074 


4 


-0.759 


1.120 


-0.122 


1.208 


-0.144 


10 


-0.785 


1.176 


-0.175 


1.310 


-0.209 



Substituting in Eq. (^) the values of ^sh determined from numerical solution of fluid- 
dynamical equations, one may calculate asymptotic parameters of the shock wave at 
given X ■ The results of this calculation are shown in Table |l|. It has been checked that 
these results agree well with jumps of flow velocity, density and pressure in Figs. 7-9. The 
last column gives the points where asymptotic velocities achieve their minima. 

4 Simulation of irradiated shocks under laboratory 
conditions 

The investigation of the effects studied in this paper in the laboratory will be possible 
by depositing energy into a target with a long deposition length. The only practical 
means for doing this is the use of heavy-ion beams. We have therefore performed realistic 
hydrodynamic simulations to estimate the magnitude of the effect with heavy-ion beams 
such as will be available with the proposed enhanced accelerator [|^ at the GSI laboratory, 
Germany. In all calculations presented in this section we consider cylindrical targets made 
of solid Gold. The targets were irradiated by singly ionized Uranium beams collinear with 
target axes. In all cases the bombarding energy of heavy ions was 500 MeV per nucleon. 

The simulations were performed using the code CAVEAT developed at Los Alamos |^ . 
The code solves the fluid-dynamical equations using the arbitrary Lagrangian-Eulerian 
[ALE) technique and the second-order Godunov-type scheme. Local values of ion energy 
deposition were calculated by decomposing the beam into a large number of beamlets and 
tracing them through the moving cells. Effects of thermal conductivity and viscosity were 
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neglected since they are not important for the considered problem. 

In the first simulation the target had 3 mm radius and 6 mm length. It was assumed 
that the beam had constant radial and temporal profiles with the beam radius = 2 
mm and the pulse duration r = 50 ns. The total penetration range of a 500 MeV/A 
Uranium ion in solid Gold of normal density is about 3.5 mm. The total deposited energy 
of the beam was 38.1 kJ. In the considered case the specific deposition rate X{x) can be 
evaluated as 

N dE 

H^) = ^-r- (39) 
vrr^r pax 

Here dE/dx is the energy loss per unit path length of a single ion in a target material 
of the density p , N is the total number of bombarding ions. The ion deposition profiles 
dE/dx have been calculated by using a realistic deposition model realized in the TRIM 
code ip. The equation of state for Gold was taken from the SESAME [|10] tables. In 
fact, we solved the fluid-dynamical equations analogous to Eqs. (0)-®, but generalized 
for the case of two-dimensional (axially symmetric) cylindrical flow. 

Two situations were compared in our simulations: irradiation with one beam from the 
right, and the same with a second beam of identical properties from the left starting with 
a 50 ns delay, i.e. exactly after the end of the first beam. The target length is chosen such 
that the second beam completely penetrates the running shock wave, and although the 
realistic deposition profile is used, the second Bragg peak is sufficiently far to the right of 
the shock. As a consequence, the local deposition from the second beam is approximately 



constant inside the shock. Fig. |T0| compares density profiles along the target axis (r = 0) 
for the two cases at several times between 50 and 200 ns. At t = 50 ns the density profile 
is, of course, not yet disturbed by the second beam. One can clearly see that the shock 
without irradiation runs towards the left, losing speed and amplitude rather slowly at t 
between 50 and 200 ns. A weaker shock caused by the backward side of the Bragg peak is 
not visible since its position is to the right of the region plotted in Fig. |1^. The decay in 
the shock strength due to the additional deposition of energy is clearly visible already at 
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t = 100 ns. The shock moves faster and the density of matter behind its front diminishes in 
accordance with the predictions of the simphfied model studied in Sect. III. At t = 200 ns 
the density jump is reduced by about 30-50%. The shock speed has increased by about 
20%, which should be clearly observable in experiments. The additional structure in the 
compressed region is apparently caused by the non-uniform deposition of the beam. 

Another possible scenario may exist even for a single beam irradiation. In the case 
when the beam has long-enough duration, the sideward expansion of the target material 
will result in depletion of density and, therefore, in a longer range of bombarding particles 
near the target axis. This enhanced range can lead to beam penetrating through the shock 
wave created at the initial stages of target irradiation. As a consequence, at late stages 
the shock wave will be weaker in a central region of the target. An example of such 



configuration is shown in Fig. |Tl|. In this case the target had a radius of 2 mm and a 
length of 5 mm. The beam had the same properties as above, except that its radial profile 
was assumed to be Gaussian with a full width at half maximum of 1.5 mm. The Gaussian 
profile was chosen not only because it is more realistic, but also because in this case the 
density starts decreasing immediately after beginning of irradiation even near the axis of 
the target. The results of calculation demonstrate how the original curved shock front 
created by the beam is destroyed in the center, leading to a ring-like structure. This 
effect can be observed after the shock wave emerges from the other end of the target. 

To check whether the effect might be observable with heavy-ion beams which will 
be available at the first stage of the upgraded GSI facility we have performed two 
additional calculations for smaller beam intensities. In the first calculation both beams 
(the first one generating the shock and the second beam irradiating it) had equal total 
energies, reduced by factor 5 as compared to previous cases. We chose the beam duration 
r = 300 ns and radius = 5 mm. The target had 3 mm radius and 6 mm length. In the 
second calculation the beam generating the shock wave was replaced by a stronger beam 
with parameters considered in the first series of simulations (see Fig. ^ and footnote on 
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page 2). In both calculations the modifications of shock velocities and compression ratios 
turned out to be at the percentage scale which is below the measurable level. Therefore, 
experimental investigation of the effects predicted in this paper requires intensive heavy- 
ion beams which will be available only at the fully upgraded GSl facility 

5 Summary and outlook 

In this paper we have investigated the behaviour of a shock wave stimulated by energy de- 
position. The case of a planar shock wave under the influence of homogeneous irradiation 
has been studied in detail. This was done by solving the one-dimensional fluid-dynamical 
equations. It is shown that the initial shock discontinuity decays into accelerated shock 
wave and a simple wave propagating into the post-shock matter. It is found that the flow 
exhibits a self-similar behaviour at late times. 

The fluid-dynamical simulations have also been performed for the more realistic con- 
ditions expected for laboratory experiments. We have studied the behaviour of a shock 
wave under the influence of additional energy deposition produced by a heavy-ion beam 
in a cylindrical target and found that the main features of the above planar solution are 
also reproduced in this more realistic situation. Another interesting effect studied is the 
weakening of the shock wave in the vicinity of the target axis due to the sideward expan- 
sion of matter and the increased deposition length of bombarding ions. These examples 
represent possible experimental set-ups for studying the effects of energy deposition on 
the structure of shock waves. Such experiments are possible at any laboratory having 
intensive heavy-ion beams with high enough deposition energies. 

In the future we are planning to study analogous processes in the case of spherical 
geometry and nonuniform density distributions that is more relevant for supernovae ex- 
plosions. Another problem is to study the influence of cooling and internal heat transport 
on the shock wave dynamics. In particular, it is interesting to investigate the role played 
by endo- and exothermic nuclear reactions during the propagation of shock front. 
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Appendix: Equations for asymptotic flow profiles 

Inserting expressions (pil|)-(|52D into Eqs. (P)-(|TT|) one obtains the coupled set of ordinary 
differential equations for asymptotic profiles V{C,), v{0 and -P(0 • 

v-^^^V' = v'V, (Al) 
v-l^)v' + {^-l)P'V = -'^, (A2) 
v-l^^P' + ilv'+l)P = l, (A3) 

where primes denote derivatives with respect to ^ . The boundary conditions for these 
equations may be derived from Eqs. (|T^-(|19]). At t ^ oo , using (|30|)-(]32D one has 

v{±oo) = 0, (A4) 
p(-oo) = V{-oo) = l, (A5) 
P(+oo) = V-i(+oo) = p2o/po. (A6) 

The above equations are trivially satisfied in the regions ^ > and ^ < ^sh • Therefore it 
is sufficient to solve (|Al|)-([7^) only in the interval (^sh, ^o) • The boundary conditions at 



^ = ^0 are given by Eqs. (|7^), (^) with the replacement +oo — > • The parameters of 
asymptotic fiow at ^ = ^sh + can be determined from Eq. (^). They are functions of 
the shock wave position . 

One can rewrite Eqs. (|A1|) -(|A^) in the vector form 

Ay' = B, (A7) 

where y is a column with the components V, v, P . The matrix A and the vector B may 
be easily obtained from the above set of equations. The derivatives y' and the matrix 
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A ^ do not exist at the points ^ where 

. (A8) 



2 



Here c = ^^oV -P^ is the scaled sound velocity. Comparing this equation with Eqs. (p!^-(p!6|) 



one can see that solutions of ( [A^ ) give the asymptotic positions of the C± and Cq char- 
acteristics. The latter is defined by Eq. (^). 

In fact, the right boundary of the nontrivial region, ^ = 5 is the C+ characteristic. 
At this point v' jumps from zero (at ^ > ^o) to some finite value = v'{^o — 0) . By 
using Eqs. (|A1|) -([A^) for small ~ ^ , one can show that their solution between the 



characteristics C+ and Cq , i.e. in the interval {^^, C,o) , is a functional of the parameters 
tty and p2o/po- At given one can solve ( |A1|) - ([A3| ) in the whole interval. Below this 
solution is marked by the subscript "+" . 

On the other hand, choosing some value of and calculating y {C,sh + 0) , one can find 
the solution of Eqs. (|Al|)-( [A3| ) between the shock front and the Cq characteristic, i.e. in 
the interval (^sh, • The corresponding solution (denoted below by the subscript "-") is 
fully determined by the parameter ■ Matching the solutions and y^ at C, = gives 
three conditions for determining three unknown parameters ^sh, ^* and a„: 

{;_(e*) = v+i^.) = l^., (A9) 
P-i^*) = P+{^*)- (AlO) 



By using Eqs. ( [A9| )-( [AT0| ) and ( |A1| )-(|A3|) it may be shown that V will be also continuous 
at ^ = ^* . 

A "shooting" method may then be used to find .^sh and asymptotic flow profiles y{^) 
at different pressure ratios x ■ These calculations agree well with the results obtained in 
Sect. IIIB. 
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Figure 1: Density profiles for different times of sliock wave irradiation. The initial pressure 
ratio X = 4 . Scales Xq and are given by Eqs. (p!2D -(|l3D. The dotted curve shows densities 
on the Cq characteristic with x(0) = 0. 
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Figure 2: The same as in Fig. |1|, but for pressure profiles. 
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Figure 3: The same as in Fig. 0, but for sound velocity profiles. 
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Figure 4: The same as in Fig. [1|, but for flow velocity profiles. 
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Figure 5: The density profiles at large irradiation times as functions of ^ = a; (At^) 
The dotted vertical lines show positions of the C± characteristics (pSf). 



-1/2 



23 




24 




Figure 7: Asymptotic profiles of scaled velocity at various ratios of initial pressure x ■ 
The thin dashed curve corresponds to the line 3^/2. The vertical dotted lines indicate 
the points of minimal velocities. Dotted lines are the same as in Fig. 
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Figure 8: Asymptotic density profiles at various x ■ The arrows indicate the values of , 
defined by Eq. (^). The dotted lines are the same as in Fig. ||. 
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Figure 10: Profiles of density along the axis of a cylindrical solid Gold target irradiated 
by a heavy-ion beam of 50 ns duration from the right at various times (the parameters of 
the beam are given in the text). The curves for times marked with an asterisk correspond 
to the situation with an additional beam of the same properties impinging on the target 
from the left after the first beam has ended. 
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Figure 11: Density plots of a Gold cylinder irradiated by a heavy-ion beam from the 
right. The darkness of shading is proportional to the density of matter. Upper and lower 
plots correspond to irradiation times 6 and 12 ns, respectively. The curved shock front is 
clearly developed after 6 ns, but is destroyed in the center already at 12 ns. The effects 
of beam penetration through the shock front are also visible at this time. Note that the 
grey scale is different in two plots: it was chosen to show geometrical effects more clearly. 
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